Shot Noise Spectrum of Open Dissipative Quantum Two-level Systems 
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We study the current noise spectrum of qubits under transport conditions in a dissipative bosonic 
environment. We combine (non-)Markovian master equations with correlation functions in Laplace- 
space to derive a noise formula for both weak and strong coupling to the bath. The coherence- induced 
reduction of noise is diminished by weak dissipation and /or a large level separation (bias). For weak 
dissipation, we demonstrate that the dephasing and relaxation rates of the two- level system can 
be extracted from noise. In the strong dissipation regime, the localisation-delocalisation transition 
becomes visible in the low-frequency noise. 
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The way a quantum two-level system (qubit) loses co- 
herence due to the coupling with a noisy environment 
has been the subject of intense research for many years 
0, 0]. This fundamental problem has received a great 
deal of attention due to recent advances in solid state 
devices in which quantum two-level systems (TLS) have 
been realized using different degrees of freedom (charge, 
spin, flux) Interest in current noise 13, in particular 
in presence of dephasing and dissipation y| , has risen ow- 
ing to the possibility of extracting valuable information 
not available in conventional dc transport experiments. 

In this Letter, we demonstrate that current noise in 
coupled quantum dots or Cooper pair boxes reveal the 
complete dissipative, internal dynamics of qubits coupled 
to external electron reservoirs. We develop a formalism 
that allows us to make quantitative predictions for the 
frequency (ui) dependent charge and current noise for ar- 
bitrary dissipative environments. We find a reduction of 
noise by coherent oscillations, weakened by increasing the 
bias or weak dissipation. The latter suppresses shot noise 
at oj = and large bias due to spontaneous boson emis- 
sion. Importantly, the dephasing and relaxation rates 
of the TLS can be extracted from noise. Our formula- 
tion includes non-Markovian memory effects [6j and the 
strong coupling limit, where we observe a re-establishing 
of the full shot noise due to the formation of polarons as 
new quasi-particles. 

In the following, we assume that the TLS is defined in 
a double quantum dot (DQD) device 0,0- We point out, 
however, that our method can also be applied to charge 
qubits realized in a Cooper pair (CP) box HOHOH 

see 

below. DQDs in the regime of strong Coulomb blockade 
can be tuned into a regime that is governed by a (pseudo) 
spin-boson (SB) model (dissipative two-level system Q), 
coupled to reservoirs [ijj H = Hsb + H res + Ht- Here, 
Hsb describes one additional 'transport' electron which 
tunnels between a left (L) and a right (R) dot with energy 
difference e and inter-dot coupling T c , and is coupled to 



a dissipative bosonic bath {TLb — Eq w Q a Q a Q)> 



'Hsb : 



■£f(» 

Q 



Q + a Q 



a z +T c a x +H. B - (1) 



The effective Hilbert space of the closed system consists 
of two states \L) = \N L + 1,N R ) and \R) = \N L ,N R + lV 
such that the system is denned by a 'pscudospin' |l2( 
a z = \L)(L\ - \R)(R\ = n L - n R and a x = \L)(R\ + 
\R)(L\ = p + p^ . The effects of the bath can be encapsu- 
lated in the spectral density J(u>) = J2q | <?q | 2 <5 — ojq), 
where loq are the frequencies of the bosons and the <?q 
denote interaction constants. When showing results we 
will be using J(u>) = 2aui [1 — Ud/u sin (uj/ujd)] e~ UJ l^ c 
for piezoacoustic phonons in lateral DQDs with aid de- 
pending on the geometry |l2j . or a generic Ohmic bath 
{tod — > 0): J(lu) = 2aLue~ l * 1 t Wc . The dimensionless pa- 
rameter a reflects the strength of dissipation and uj c is a 
high energy cutoff • The coupling to reservoirs H re s = 



Cfe 



is described by Ht — J2k ^Vk^k 



H.c ), where s a = |0)(a| (a=L,R) and the extra state 
|0) = \Nl,N r ) describes an 'empty' DQD, such that 
1 = no + n L + n R . 

The full model described by H allows to study non- 
equilibrium properties, such as the inelastic stationary 
current or current noise, through an open dissipative 
TLS. We describe its dynamics by a reduced, with respect 
to reservoirs, statistical operator p(t). Introducing the 
vectors A = (n L , h R ,p,p^) T , T = (r L ,0,0,0) T = r L ex 
and a matrix memory kernel M, the equations of mo- 
tion (EOM) of the expectation values [13| (with (O) = 



E 



i=0, L,R 



Trbath(*|Op(£)|i)) read in matrix form 



(A(i)) = (A(0))+ / dt'{M(t-t')(A(t'))+r}. (2) 



Eq. (0) can be solved in Laplace space as (A(z)) = \z — 
zM(z)] _1 ((A(0)) + T/z) and serves as a starting point 
for the analysis of stationary (1/z coefficient in Laurent 
series for z — > 0) and non-stationary quantities. The 
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memory kernel has a block structure 

zM(z) = 



-G 



G = 







(3) 



where T c = — iT c (l — a x ), and the coupling to the reser- 
voirs within Born and Markov (BM) approximation with 



respect to Ti T [13,113 is S iven b Y r « = ^ J2k a \ V k \ 5 ( e ~ 
e k a ) (we assume Fermi distributions for the reservoirs 
/l = 1 and fn = 0; large voltage regime). The blocks 
D z and S z are determined by the EOM for the coherences 
(off-diagonal elements) (p) = (?r )* and contain the com- 
plete information on dephasing of the system. In general, 
no exact solution is available but we will present approx- 
imate results now: for weak coupling to the bosons, one 
can use perturbation theory (PER) in a in the correct 
basis of the hybridised states of the TLS. In BM approx- 
imation, the resulting expressions are: 



D 



PER 



7+ -7- 
7+ -7- 



^PER 



E 
E* 



(4) 



where E = ie - j p - j p = 2tt^ J (A) coth (/3A/2) 
and 7± = J(A) coth (/3A/2) =p J(A) com- 

pletely determine dephasing and relaxation in the sys- 
tem. Here, A = e 1 + 4T C 2 is the hybridization splitting 
and @ = l/k B T. 

On the other hand, for strong electron-boson coupling, 
one has to start from a polaron-transformed frame 
(strong coupling, POL), leading to an integral equation 
|l2| which involves the boson correlation function C (t) = 



J(u) 



(i 



Introducing C e (z) : 
resulting matrices in z-space are 



cosw*) coth(^ 
= J °° die 



D 



POL 



lT r 



-I C*_JC e 

1 -C- e /Cl 



'-.POL 



+ i sin u)t 



E 
E* 



the 



(5) 



with £[*] = z - l/d* ] (z) - T R /2. In contrast to the 
PER solution, where M(t) = M = zM(z) is time- 
independent, M pol (t) is time-dependent and zM{z) de- 
pends on z in the POL approach ^ij. We note that 
Re[C e (z)]| z= ± i(i) = nP(e =F where P{e) is the proba- 
bility for inelastic tunneling with energy transfer e Q- 

As mentioned above, our model describes a CP box 
as well, the transport through the DQD being analo- 
gous to the Josephson Quasiparticle Cycle (JQP) of the 
superconducting single electron transistor (SSET) with 
Ec S> Ej, such that only two charge states, |2) (one ex- 
cess CP in the SSET) and |0) (no extra CP), are allowed. 
Two consecutive quasiparticle events (with rates T2 and 
Ti) couple |2) and |0) with another state |1) through the 
cycle |2) — * |1) — * |0) |2). Interdot tunneling is anal- 
ogous to coherent tunneling of a CP through one of the 
junctions, and tunneling to and from the DQD is anal- 
ogous to the two quasiparticle events through the probe 
junction in the SSET \Ti\ . 



Current noise, which is described by the power 
spectral density S/(oj) = 2 dre' 1 " 7 <Si(r) = 

E» dre^({Ai(r), AJ(0)}>, with A/(*) = /(*) - (/(*)), 
is a sensitive tool to study correlations between carriers 
4J. The Fano factor (7 = S J,^ ) quantifies deviations 
from the Poissonian noise, 5/(0) = 2ql, which charac- 
terizes uncorrelated carriers with charge q. Importantly, 
Si(uj) has to be calculated from the autocorrelations of 
the total current I(t), i.e. particle plus displacement cur- 
rent Q. Using current conservation together with the 
Ramo-Shockley theorem, I(t) = aljt) + blu(t) (a and b, 
with a + b = 1, depend on each junction capacitance |j|), 
one can express iSr(w) in terms of the spectra of particle 
currents and the charge noise spectrum Sq(u>) |16|. 

Si(cj) = aSi L (uj) + bSi R (uj) - abLu 2 S Q (uj). (6) 

Note that in symmetric configurations, a ~ &, the charge 
noise reduces the contribution from particle currents to 
the noise spectrum. For a ss 1 or b ss 1 the main con- 
tribution to noise comes from particle currents. At zero 
frequency 5/(0) = Si L (0) = Si R (0). Sq(uj) is defined as 



S q {lu) = lim / dTe WT ({Q(t),Q(t + T)}} 

J —OO 

= 2Re {f(z = iu) + f(z = -iu)}, (7) 
where Q = ul +nn and f(z) is the Laplace transform of 



f(r) 



E < 

i,j=L,R 



hi(t)hj(t + t)) 



(8) 



and can be evaluated with the help of the charge cor- 
relation functions C Q (r) = (n a (t)A(t + r)), as /(r) = 
(e 1 +e 2 )[C L (T)+C jR ,(r)]. The EOM for C Q (r) can be ob- 
tained from the quantum regression theorem |l7| whose 
solution is again expressed with the help of the resolvent 
[z-zM(z)}- 1 , cf. Eq. ©. 

To calculate the contribution of particle currents 
to noise, we need to relate the reduced dynamics of 
the qubit described by Eqs. (2-3) to reservoir op- 
erators. For Si r (lu), we introduce the number n of 
electrons that have tunneled through the right barrier 
lis! flf^ which defines generalized expectation val- 
ues as = I]i=o.L.fl Tr bath("-, i\6p(t)\n, i) (such that 

(6) = E„° (n) ) and write 



■ («) 
n L/R 



Rn 



(n-l) 



±^L/Rn { Q ] ± iT c 



,(«) _ 



)lt 



(9) 



and correspondingly for and [pW]t and the left bar- 
rier. Eqs.® allow one to calculate the particle cur- 
rent and the noise spectrum from P n (t) — Uq (t) + 
nf^ (t) + (t) which gives the total probability of find- 
ing n electrons in the collector by time t. In particular, 
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FIG. 1: a) Fano factor vs. bias e for different dissipative 
couplings a. Parameters T c = 3, V — 0.15, u) c = 500, 
u>d — 10, T = 2 (in iieV) correspond to typical experi- 
mental values in double quantum dots. Lines: acoustic 
phonons, circles: generic ohmic environment ujd — (see 
text), b) Frequency dependent current noise (a = 0, T = 0, 
r = 0.01). Inset: (Top) Contribution to noise from parti- 
cle currents Si R (ui)/2eI. (Bottom) Charge noise contribution 
ui 2 S Q (Lo)/8eI. a = b = 1/2. 



Iri(t) — e J2 n n Pn(t) and Sj R can be calculated from [2(j 
S Ir (uj) = 2uje 2 I dtsin^t)- [(n 2 (t)) - (t(I)) 2 } ,(10) 



where j- t {n\t)) = £„n 2 P„(i) = TrYZ*™^® 



n=0 11 R 



(n) 



(t). Solving Eqs.JJjJ with the initial condi- 
tion n R (0) = S n .on R (0), where n R (0) is the stationary 
solution of Eqs. (2) yj|, we get 



Si r (lj) = 2eI{l + T R [n R (-iuj) + n R {iu)]}, (11) 

with zh R (z) = T L g + (z)/N(z), where N(z) = [z + T R + 
g-(z)](z + T L ) + {z + T R + T L )g+{z) and 



9+[-\{ z ) = ±iT c (ei - e 2 ) 



D z e x \ 



(12) 



Eqs. I|11I12I) demonstrate the dependence of the current 
noise on the dephasing via the two-by-two blocks D z and 
t z , cf. Eq. Explicitly, 

T c ( lp + T R /2 + z)- e 1± 



gl ER (z) ee 2T/ C ' '' 



,POL 



Tz 



( 7p + r fl /2 + zf + e 2 
cW ]£ (z) 



+ (C C*) 



(13) 



A similar derivation yields S 1 ^ (w) = (w) . The explicit 
expressions Eq. (|11I12I13I) , together with the inverse of a 
4 by 4 matrix for the charge noise Eqs. Q, yield our key 
quantity Si(u>), Eq. 

Zero frequency (Shot Noise).- In the zero frequency 
limit z — > 0, one obtains 



5j(0) = 2el ( 1 + 2r fl — [zn R (z)] z=0 



(14) 




FIG. 2: Effect of ohmic dissipation on current noise near reso- 
nance (e = 10, T = 0.01, and a = 0.005, 0.01, 0.02 correspond- 
ing to 7 P ss 4. 74r, 9. 47r, 18.95r). Inset: (Right) pseudospin 
correlation function S z (uj). Arrows indicate relaxation rate 
( r + 7p)/ a ~ 0.005 for a = 0.005. (Left) low frequencies 
region near shot noise limit w = 0. 



Eqs. (|14|) allows to investigate the shot noise of open 
dissipative TLS's for arbitrary environments. In con- 
trast to non-interacting mesoscopic conductors, the noise 
cannot be written in the Khlus-Lesovik form 5/(0) = 
2e 2 J j^-t(E)[l — t(E)] with an effective transmission co- 
efficient t(E). Without bath, we recover the results of 
Ref. (shot noise of DQD's) and Ref. [H| (shot noise of 
the CP box). For a = and T = T L = T R (Fig. la, solid 
line), the smallest Fano factor is reached for e = where 
quantum coherence strongly suppresses noise. The max- 
imum suppression (7 = 1/5) is reached for T = 2y/2T c . 
For large e > (e < 0) the charge becomes localized 
in the right (left) level, 5*7(0) is dominated by only one 
Poisson process, namely the noise of the right (left) bar- 
rier, and 7 — > 1. For a ^ 0, this mechanism is strongly 
affected by the possibility of exchanging energy quanta 
with the bath. The effect of a bosonic bath on noise in a 
mesoscopic scatterer was first discussed in [il]. For the 
TLS discussed here, spontaneous emission (for e > 0) oc- 
curs even at very low temperatures 0, an d the noise 
is reduced [j| well below the Poisson limit (Fig. la). The 
maximum suppression is now reached when the elastic 
and inelastic rates coincide, i.e., 7p = T R , as we have 
checked numerically. For large couplings, spontaneous 
emission leads to a very asymmetric Fano factor that 
goes from 7* 1 to 7 « 0.5 ase changes sign (not shown 
here). 

Finite frequencies - For finite w, the numerical results 
for a — are shown in Fig. lb where we plot Si(u>) 
for different values of e. The background noise is half 
the Poisson value as one expects for a symmetric struc- 
ture. 7 deviates from this value around uj = where 
the noise has a peak and u) = A where the noise is sup- 
pressed. The dip in the Fano factor directly reflects the 
resonance of the subtracted charge noise Sq(uj) around 
A (inset Fig. lb), cf. Eq. JJJJ. An increase of e 
localizes the qubit and, thus, the zero-frequency noise 
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FIG. 3: Low frequency current noise in qubit with strong 
ohmic dissipation (T = 0, e = 10, T c = 3, Fl = = 0.01). 
Inset: Shot noise for u — 0. 



I(|e|) oc ImC_| £ | through a strongly dissipative quantum 
ring containing a quantum dot with bias e 22]. Although 
POL becomes less reliable for a < 1 and smaller bias, the 
non-symmetry in e of the shot noise and the inelastic cur- 
rent oc ReC e reflects the 'open' topology of our TLS in 
the non-linear transport regime. 

To conclude, our results demonstrate that frequency- 
dependent current noise provides detailed information 
about the internal, dissipative dynamics of open quan- 
tum two-level system such as double quantum dots or 
Cooper pair boxes. The weak coupling regime should be 
close to current experiments 23] in these systems, where 
we expect our predictions to be tested in the near future. 

We thank Markus Biittiker, Leo P. Kouwenhoven and 
Till Vorrath for useful discussions. This work was sup- 
ported by EPSRC GR/R44690, DFG BR 1528 and by the 
MCYT of Spain through the "Ramon y Cajal" program 
and grant MAT2002-02465 (R. A.). 



reaches 7 — ► 1. Moreover, the dip in the high frequency 
noise at u> = A (Fig. lb) is progressively destroyed (re- 
duction of quantum coherence) as e increases which is 
consistent with the previous argument. A similar re- 
duction of the dip at u> = A occurs at fixed e and T 
with increasing dissipation (Fig 2) in the weak coupling 
(PER) regime. This behavior demonstrates that Si(ui) 
reveals the complete internal dissipative dynamics of the 
TLS. The above argument can be further substantiated 
by plotting also the symmetrized pseudospin correlation 
function S z (u) = 1/2 due*" (R(r), a z }) (Fig. 2, 
right inset) which is commonly used to investigate the 
dynamics of the SB problem |2j. Both functions reflect 
in the same fashion how the coherent dynamics of the 
system progressively gets damped by the bosonic bath. 

In particular, the dephasing rate can be extracted from 
the half- width of Si(u) around u = A. For an Ohmic 
environment, 7^ = 7 P /2 + 2Tra(-^) 2 kBT, such that the 
total dephasing rate is j d (T = 0) = 7^ + r/2 = (7 P + r)/2 
(Fig. 2, arrows denote full-width, i.e. 27^ w j p as a 
increases). Close to u = 0, the peak in Si(u) for a = 
changes into a dip around u = reflecting incoherent 
relaxation dynamics for ci / 0. The half-width is now 
given by the relaxation rate such that the full-width of 
Si(u>) around u> — is twice that of the high frequency 
noise (Fig 2, left inset). 

The results for the strong coupling (POL) regime are 
presented in Fig 3. Near u = 0, POL and PER yield 
nearly identical results for the noise Si(u>) at very small 
a (not shown here). The cross-over to Poissonian noise 
near u> = with increasing a indicates the formation of 
localized polarons. The delocalisation- localisation tran- 
sition 0, 13] of the spin-boson model at a = 1 is reflected 
in a change of the analyticity of C £ and the shot noise 
near zero bias (Fig 3, inset). Similar physics has been 
found recently in the suppression of the persistent current 



[1] A. J. Leggett, et al , Rev. Mod. Phys. 59(1), 1 (1987). 

[2] U. Weiss, Quantum Dissipative Systems, (World Scien- 
tific, Singapore, 1993). 

[3] Y. Nakamura et al, Nature, 398 (1999), 786; D. Vion 
et al, Science 296, 886 (2002); I. Chiorescu et al, ibid. 
299, 1869 (2003); Yu. A. Pashkin et al, Nature 421, 823 
(2003). 

[4] Y. M. Blanter and M. Biittiker, Phys. Rep. 336, 1 (2000). 
[5] A. Shimizu and M. Ueda, Phys. Rev. Lett. 69, 1403 

(1992) ; Phys. Rev. Lett. 8 6 3694 (2001). 

[6] D. Loss, D. P. DiVincenzo, cond-mat/0304118 (2003). 

[7] T. Hayashi et al., Phys. Rev. Lett., 91, 226804 (2003). 

[8] T. Fujisawa et al , Science 282, 932 (1998). 

[9] Yu. Makhlin et al , Rev. Mod. Phys. 73, 357 (2001). 
[10] A. A. Clerk et al , Phys. Rev. Lett. 89, 176804 (2002). 
[11] M.-S. Choi et al, Phys. Rev. B 67, 045105 (2003). 
[12] T. Brandes and B. Kramer, Phys. Rev. Lett. 83, 3021 
(1999). 

[13] T. Brandes and T. Vorrath, Phys. Rev. B 66, 075341 

(2002) . 

[14] For the SB problem, POL is equivalent to the double- 
path integral non-interacting blip approximation 0, 
and works well for e = and e ^> T c . 

[15] T. H. Stoof and Yu. V. Nazarov, Phys. Rev. B 53, 1050 
(1996). 

[16] D. Mozyrsky et al, Phys. Rev. B 66, 161313 (2002). 
[17] See, e.g., H. J. Carmichael, An Open System Approach 

to Quantum Optics, (Springer- Verlag, Berlin 1993). 
[18] B. Elattari and S. A. Gurvitz, Physics Letters A 292, 

289 (2002). 

[19] R. Ruskov and A. N. Korotkov, Phys. Rev. B 67, 075303 

(2003) . 

[20] D. K. C. MacDonald, Rep. Progr. Phys. 12, 56 (1948). 
[21] M. Ueda and A. Shimizu, J. Phys. Soc. Jap. 62, 2994 

(1993) . 

[22] P. Cedraschi and M. Biittiker, Annals of Physics 289, 1 
(2001). 

[23] R. Deblock et al , Science 301, 203 (2003). 



